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Abstract 

It is shown that theories already presented as rigorous mathematical formalizations of 
widespread manipulations of Dirac's delta function are all unsatisfactory, and a new 
alternative is proposed. 
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I. Introduction 



There is a "myth" in Mathematical Physics according to which: 

The widespread manipulations of Dirac 's delta function can be rigorously formalized 
by the Theory of Distributions. 

It is easy to verify, however, that this statement can be taken as true only by someone 
who does not know the nature of such manipulations or by someone who does not know 
the available theories about the concept of distribution. 

The most elementary manipulations of the delta function, and also the most often 
found in literature, involve the composition of S and a real function. For example: 



where ai, 02, . . . , fin are the simple roots of a real function g . 

Considering the delta function as a distribution according to the best known definition 
of this concept ^ — linear functionals on some test functions space — the above formulae 
could not be written as such, since the composite of a distribution and a generic real 
function simply cannot be done. 

Using the sequential approach ^ we can define the composition of a distribution and 
some very simple real functions, but the function which sends x to (x"^ — a^) is not among 
them. Therefore, also in this context, not even the first formula above could have been 
written. 

It is clear that we can use the above equality to "define" the symbol ^[(/(x)] , but 
in any way could this be considered a "formalization" of the relatively simple idea which 
one has in mind when writing it down. It is obvious that the formal meaning ascribed to 
the basic and elementary manipulations of Delta Calculus must be related to the mental 
process usually understood in these manipulations. 

In this sense, the sequential approaches to the concept of distribution are much su- 
perior to that of Schwartz-Sobolev, which systematically uses the "theorems" of Delta 
Calculus as "definitions" of the involved concepts. For instance: the distributions are de- 
fined as linear functionals on some test functions space so we can characterize the delta 




(a > 0), 



or, more generally: 
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function by its sifting property: 



/ S{x)f{x)dx = f{0) 



J — oo 



That procedure, considered "elegant" by many mathematicians, merely dismisses the fact 
that the sifting property itself is a basic result of the Delta Calculus to he formally proved. 

Dirac has used a simple argument, based on the integration by parts formula, to get 
the sifting property of the derivative 5' of the delta function: 



The theory of distributions as linear functionals, instead of defining the integral of a dis- 
tribution and so proving that it satisfies some kind of integration by parts formula, just 
uses the formula deduced by Dirac for the delta function: 



as "distributions derivative definition" . 

Proceeding systematically in this fashion, it is not surprising that some "unnecessary" 
constructions must be ignored after certain symbols and operations have been "defined" , 
because we know Delta Calculus has a nontrivial mathematical content. 

The real surprise is the consideration of this Theory of Distributions as a satisfactory 
formalization of Delta Calculus, despite its ignoring the basic operational rules of this 
Calculus. 

All proposed formalizations of Delta Calculus fail because they try to "generalize" 
the real functions without having generalized its variables. This will become clear if we 

compare those formalizations to the one proposed in this work. 

In Ref. 4 we presented an extension process that can be applied to any set. This 
process was applied ^ to the ordered field of real numbers, formally defining what we call 
virtual numbers. Among those numbers we have infinitesimal and infinite quantities, which 
were so used in an reorganization of Differential and Integral Calculus. The same extension 
process was applied ^ to the set of real functions, introducing what we call virtual functions. 
So, the basic techniques of Infinitesimal Calculus were generalized to these mathematical 
objects. 

The concept of "virtual function" is very close to the original idea of "improper func- 
tion" created by Dirac: the limit of a sequence of real functions ^ . The essential difference 
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is that Dirac did not have access to "improper numbers" which could have been used as 
argument of these functions, whereas the virtual functions can be evaluated at (previously 
defined) virtual numbers, exactly as we do with real functions and numbers. 

We will here use many concepts and results presented in Refs. 4, 5 and 6, including 
the notational conventions adopted then (like the notation for reduced integrals stated in 
Ref. 6). 

This paper begins by defining, in Sec. II, the class of virtual functions having the 
characteristics supposed by Dirac for his delta function. In Sec. Ill we prove the "Dirac's 
functions" sifting property. Section IV is dedicated to discussing the meaning of the 
operational rules which make up Delta Calculus. In Sec. V we prove the basic rule for 
compositions involving delta function (above). Other well known results of Delta Calculus 
are demonstrated in the last two sections. 

II. Dirac's Functions 

The aim of this section is to characterize the virtual functions which possess the basic 
properties assumed by Dirac for his "delta function" , and to show that there exist virtual 
functions in .F(R) CJ^(R) with those properties. 

We will say that a virtual function S e .F(R) is a Dirac 's function when: 

(i) 5 is defined at any ^ G R and is non-negative: 

5(e) >0, for every R ; 

(ii) 6 is integrable with 




S{x) dx = 1; and 



(iii) there exists a positive infinitesimal £ e R such that 

|e|>e ^ 5(0 = 0. 

Before proceeding, it is convenient to verify that there exist virtual functions under 
those conditions. This is not so difficult: we know there exists a real function /: R — > R 
infinitely derivable which vanishes outside the interval (—1,1), increasing between — 1 
and 0, where it takes the value /(O) = 1 , decreasing between e 1, and with unitary 
integral: 

J f{x) dx = 1. 
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From such a function / , we can define a virtual function ^q: R — > R by: 

(5o(C) = oo/(ooO, 



i.e., So e .^(R) is the virtual function represented by the sequence (/i, /2, . . .) G S[J^(R)] 
given by: 

fn{x) = nf{nx). 

It is easy to see that this virtual function satisfies the three conditions above: 

^o(0 for every ^ G R , 

So{x) dx = 1, 



/ 



and 

\^\>d ^ 5o{0 = 0- 

Therefore 60 is a Dirac's function. 

That virtual function do perfectly corresponds to the image given by Dirac ^ : 
"To get a picture of S{x) , take a function of the real variable x which vanishes 
everywhere except inside a small domain, of length e say, surrounding the origin x — , 
and which is so large inside this domain that its integral over this domain is unity. The 
exact shape of the function inside this domain does not matter, provided there are no 
unnecessarily wild variations (for example provided the function is always of order ). 
Then in the limit £ — > this function will go over into 6{x) ." 

Moreover, the function (5o : R — > R constructed above is infinitely derivable (so 
continuous), normalized in the strong sense: 



<^o(o^e = i, 

and, evaluated at real values of its argument, provides: 

, / V ( 00 , if a; = 0: 
'^°(^) = l0 ,ifx^O. 

In order for a virtual function 5: R — > R to be a Dirac's function, it is not necessary 
that 5{0) = 00 , or even that 5{0) 00 . But the condition 5{x) = , for every real x ^ , 
immediately follows from the existence of an infinitesimal s such that |^| > e 5{^) = . 
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For instance, the virtual function 7/;: R — > R given by: 

00 



1 + Oo2^2 

is such that ipi^) ~ ^oi every x ^ : 

= ^— ^ — ^ ~ 0. 
o + oox^ 

But V'lO 7^ 5 C ^ 7 ^-iid so -0 is not a Dirac's function. (That function is the 

derivative of = arctan(oo^) , and it was used in various examples in Ref. 5.) 



It should be clear that there exist many distinct Dirac's functions in J^(R) . For 
example, if 5o : R — ^ R is the function above constructed then the virtual functions 
: R — > R and 5_ : R — > R given by: 

and 

are two other infinitely differentiable Dirac's functions distinct from Sq . The function S+ 
is the class of the sequence {gi,g2, • • •) ^ ^[-^(I^)] defined by: 

2 

9n{x) = fn{x ), 

n 



where fn are the above functions which represent Sq e ^(R) • Analogously, 5_ is the 
class of the sequence (/ii, /i2, . . .) ^ ^[-^(R)] defined by: 

2 

hn{x) = fn{x+ -). 

n 

A Dirac's function like S+ can be convenient if we intend, for instance, to apply the 
Laplace transformation, since it vanishes for every ^ < , but is normalized in the sense: 

/•oo 

/ S+{x) dx = 1. 
Jd 

Furthermore, it is easy to verify that: 

If Si and S2 are two Dirac's functions then the virtual function ^3: R ^ R defined 

by: 
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is also a Dirac 's function. 

We have not required a Dirac's function to be always continuous or derivable, so we 



can take "square pulses" like the function 5n £ -^(R) represented by the sequence: 

/ \ _ / "'/2 , if \x\ < 1/n; 
gn{x)-<^^ ,if|x|>l/n. 

This Dirac's function, which is discontinuous at ^ = —d and ^ = d , can also be useful in 
specific situations. It is such that: 



00/2 ,if|e|<a; 

,if|el>5- 



III. Sifting Property 



Our aim in this section is to demonstrate the sifting property of Dirac's functions. We 
will say that a real function is continuous around a; e R when there exists a real open 
interval containing x in which it is defined and continuous. 

If S: R ^ R is a Dirac's functions, and /: R — > R is a real function continuous 
around the origin, then: 

poo 

S{x)f{x)dx = f{0). 



' —00 

Proof: First, we note that there exists a positive infinitesimal e such that: 

\^\>e ^ 6(0 = 0, 

for 5 is a Dirac's function. By the additivity with respect to the virtual integration 

interval ^ , we have: 

/ — e /■£ poo 

-00 J —s Js 

On the right-hand side of the above equation, the integrands of the first and third integrals 
vanish, for any function / . So those integrals exist and are equal to zero. In the second 
integral (right-hand side) we have the product of two integrable functions, so the left-hand 
side integral exists and: 



/oo re 
mfiOd^= / mfiOd^ 
-00 J —s 



It is clear that, in R , the following statement holds: if g and h are two real functions 
defined between —a and a , where a is a positive real number, with g integrable and 
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non-negative and h continuous, then there exist two real numbers ci and C2 between 
—a and a such that: 

/a pa pa 

g{x) dx < / g{x)h{x)dx < h{c2) / g{x)dx. 
-a J — a J — a 

So we have, by the Virtual Extension Theorem (VET, Ref. 4), that if and ip are two 
virtual functions defined between —a and a , where a is a positive virtual number, with 
integrable and non- negative and V continuous, then there exist two virtual numbers 
7i and 72 between —a and a such that: 

Hii) r < r mmd^ < C md^. 

J —a J — a J —a. 

Now making a = e , (f) = S and i/j = f , we conclude that there exist two virtual numbers 
7i and 72 between —e and s such that: 

/(7i) £ dc < f mfio dc < /(72) £ dc 

Since 7i and 72 are between the infinitesimals —e and e , they are also infinitesimals 
themselves. Then it follows from the continuity of / at the origin that: 

fill) ^ /(O) ^ /(72), 

and therefore: 

/oo 
mfiOd^^fio). ■ 
-00 

It is important to note that the sifting property demonstrated above holds for all 
Dirac 's functions. If we are working with a particular Dirac's function in a specific context 
then we can deduce variants of this property by analogous methods. For instance, we 
can use the function from the previous section to "unilaterally" sift a real function / 
defined only for positive values of x . 

IV. Dirac's Equivalence 

Delta Calculus is based on a set of operational rules, whose meaning has been estab- 
lished clearly by Dirac himself^ : 

"There are a number of elementary equations which one can write down about S func- 
tions. These equations are essentially rules of manipulation for algebraic work involving 
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S functions. The meaning of any of these equations is that its two sides give equivalent 
results as factors in an integrand." 

We will say that two integrable virtual functions 0: R — > R and 7/;: R — > R are 
equivalent when, for every continuous real function /: R — > R , we have: 



It is important to note that, according to this definition, both integrals must be reducible 
in order for the two functions to be equivalent. That means, if 0: R — > R is such that 
there exists a continuous real function /: R — > R which makes the virtual integral 



irreducible, then is equivalent only to itself. For example, if 5 is a Dirac's function and 
A; e R a non-zero real constant, then the sum function (S + k) is equivalent only to itself, 
since the integral 



is not reducible. 

We will indicate that two virtual functions (f) and are equivalent by writing: 



for every Dirac's function 6 (this is a simple consequence from the "composition rule" 
which will be stated and proved in the next section). We use a real variable in this notation 
(Latin, not Greek letter) to reinforce that the two expressions are interchangeable only in 
reduced integrals. 

It is not difficult to verify that " = " is an equivalence relation on the set of integrable 
virtual functions defined at any ^ e R . This relation will be called Dirac's equivalence. 

The above example can be used to gain a fiash of intuition about the nature of this 
equivalence: the graph of function 5{2x) is "the same height as one of the function 5 , 
but half of its width"; whereas the graph of function {l/2)5{x) is "the same width as the 






0(^) 



ip{x). 



For instance: 
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one of function S , but half of its height". Thus, although fairly different, those functions 
are equivalent as factors in an integrand: 



since both delimit "the same area infinitely concentrated around the origin" , whereas any 
continuous real function / "varies slowly" . 

To discuss the equivalence of Dirac's functions, let us first consider the following 
definition: an integrable virtual function </>: R — > R is sifting when: 



for any continuous real function /: R — > R . 

It is not difficult to see that the set of all sifting functions is an equivalence class 
according to the Dirac's relation. 

Every Dirac's function is sifting (as we saw in the previous section), but not every 
sifting virtual function is a Dirac's function. As a counterexample, we can take the sequence 
ifi: /2, • ■ ■) ^ ^[-^(R-)] ^sed in Sec. I to define Sq , and alter all fn at one and only one 
point: 



This sequence {gi,g2,...) G S[.F(R)] represents an integrable virtual function 0: R — > R 
which is clearly sifting, but not a Dirac's function, for ■^(7) = 3 . 
Even so, we can affirm: 

Any two Dirac 's functions are always equivalent to each other. 

This fact makes two conventions which are part of the traditional language of Delta 
Calculus compatible: 

(i) All Dirac's functions are represented by the same symbol: " 5 " . 

(ii) The Dirac's equivalence is indicated simply by the equality symbol: " = " . 
Those conventions drastically simplify notation, and certainly do not jeopardize the 

rigour of a scientific work if used properly. Nevertheless, since they might generate some 
confusion, we will continue to explicitly distinguish equality ( = ) from equivalence ( = ), 
and to represent distinct Dirac's functions by distinct symbols. 

The traditional language of Delta Calculus also does not establish a clear distinction 
between sifting functions and Dirac's functions, since the latter are generally handled using 
the sifting property. That is another possible source of confusion, which we should keep 
in mind, mainly when the above conventions are used. 





oo 




fn{x) , ff a; 7^ 7; 
3 , if a; = 7. 
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V. Composition Rule 

Dirac's functions, as defined in Sec. II, are virtual functions 6 G -^(R-) with certain 
specific characteristics. Therefore, it is clear that they can be composed with any other 
virtual function (f) . Besides, the composite virtual function {So(f)) is defined at every point 
in the domain of (p i since Dirac's functions are defined at any virtual number ^ e R . 

Several of the operational rules of Delta Calculus involve the composition of a Dirac's 
function 5 and a real function g . We looked at the meaning of those rules in the previous 
section: they establish Dirac's equivalences for the composite virtual function {5 o g) . To 
discuss these equivalences, we will suppose that the real function g is defined on the whole 
virtual extension R of the real line. 

Dirac's functions do not vanish only around the origin, so: 

7/ 5: R ^ R is a Dirac's function and ^: R — > R a real function for which there 
exists a positive r e R with \g{x)\ > r for every x & R , then the composite function 
(5 o gf): R — > R is identically null: 

^[giO] = 0> for every ^eK . 

In this case, it is cleat that ^[^'(O] = • For instance: 

5{x^ + l) = and 5{smx + 2) = 0. 

Let us consider now a real function whose image approximates the origin: g{x) = x^ , 
for example. If we take the Dirac's function 5- defined in Sec. II: 

5_(e) = 5o(e+2a), 

then: 

5_(e') = 0, for every R, 
since d-{^) — for any ^ > . So it is clear that: 

S_{x'^) = 0. 

On the other hand, composing the "square" Dirac's function 5n (also defined in 
Sec. II) and that same function g{x) = , we get a non- vanishing composite: 

\^\<Vd^Sn{e) = Y- 
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Furthermore, it easy to see that: 
for: 

/oo /"VS 
-oo V-Va ^ ^ 

So, we have: 

which shows that the equivalence class of this composite depends on the particular Dirac's 
function chosen. That means, there is not a virtual function such that 

for any Dirac's function 5 . Therefore, there is not an operational rule for the compos- 
ite 5{x'^) in the traditional language of Delta Calculus, which deals only with generic 
equivalences that hold for all Dirac's functions. 

This dependence of the class of the composite with respect to the selection of the 
Dirac's function might also occur if the graph of the real function ^: R — > R is asymptotic 
of the x -axis. As an example, for g{x) — we have: 

5-{e^) = for every ^ e R , 

whereas: 

S+{e^) ^ 0. 

However, the class of the composite function will not depend on the chosen Dirac's 
function if the graph of the real function ^: R — > R only "crosses" the a; -axis in a finite 
number of roots, not approximating it in any other region: 

Let 5: R — > R be any Dirac's function and ^r: R — > R a real function with a finite 
number of roots ai, a2, . . . , On i o-i^d such that there exist n+1 positive real numbers 
■ • • jT^n having the following properties: 

(i) the n intervals [aj — ri,ai + r^] are disjoint in pairs and 

\9{x)\ > r, 

for any a; e R outside those intervals; and 
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(ii) inside the n intervals [oj — r^, + rj] , the function g is differentiable and its 
derivative does not vanish. 

Then, for any real function /: R — > R continuous on the n intervals [a^ — r^, + rj] , 
we have: 

Proof: Since the composition vanishes outside the n intervals [a^ — r, + r] , we 

have: 

-oo J tti—ri 

We will calculate the integrals in the above sum separately. First, we note that, for each 
i = l,2,...,n, the restriction of the function g to the interval [oj — r^, Oj + r^] admits an 
inverse hi : 

hi[g{x)] = X, for every x e [ai - rj, Oi + rj] , 

and this inverse hi is monotonic, differentiable and its derivative does not vanish between 
gi^i + ri) and g{ai - r^) . 
Thus, changing variables: 

I^ = 9i0i ^ = hi{ii) and dii = g'{^)d^, 



we get: 

If g'{ai) > then g{ai — ri) < < 5'(aj + ri) , so the sifting property guarantees that: 

On the other hand, if g'{ai) < then g{ai + rj) < < g{ai — ri) , so the sifting property 
provides: 

Thus, in any of both cases we have: 



/ smmod^ 

J ai—Ti 



fjaj) 
\9'{ai)\' 
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For any number a e R , the function g{x) = x — a satisfies the hypothesis of the 
above result, so we have the well known sifting property of translated Dirac's functions: 



/oo 
S{x - a)f{x)dx = f{a), 
-oo 



for any /: R ^ R continuous around a G R . 

So, we conclude that, for any real function g under the above conditions, and for any 
Dirac's function S : 
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This formula will be called composition rule. 

The sifting property of translated Dirac's functions implies: 

f{x)5{x -a) = f{a)5{x - a), 

for any function /: R — > R continuous around a e R . 
The composition rule, on its turn, shows that: 

S{ax) = ^^S{x) (a 0), 



and 



S{x'^ -a"^) = ^ [S{x -a) + d{x + a)] (a > 0), 



for any Dirac's function 5 . 

VI. Contraction of Dirac's Functions 

Dirac's functions, as defined in Sec. II, are integrable virtual functions 5: R — > R . 
Thus, it is clear that, for any /? e R , the virtual function which assigns — (3) to ^ is 
also integrable. The result below will be used as a lemma in the following demonstration: 

7/ 5: R — > R is a Dirac 's function, and e R a finite virtual number, then: 



/oo nOO 
5{^-p)d^= / 5(0 d^, 
-oo J — oo 



Proof: Changing the variables n = ^ — (5 we get: 



/_ 



>3 l>00 — /3 

oo J— oo— 13 

— oo foo poo—/3 



/oo i>oo—p 
5{ii)diJi-\- i S{iJ,)diJ,, 
-oo J oo 



00-/3 
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but the integrands of the above first and third integrals vanish, for P is finite and 5 is a 
Dirac's function. I 

If Si and S2 are two continuous Dirac's functions then the virtual function ^3: R — > 
R defined by: 

j>00 

m)- / 5i(e-/9)52(/9)rf/? 



is also a Dirac 's function. Besides, we have: 

/oo 
-00 

for any finite a e R . 

Proof: It is quite easy to see that 53(^) > , for every ^ e R , since the above integrand 
is non-negative. To show that ^3(0 vanishes far from the origin, we first note that there 
exist positive infinitesimals £1 and £2 such that: 



iei>£i^5i(0 = o 



Taking £3 = £1 + £2 we get: 



and 



1^1 >£2^ 52(0 = 0. 



|e|>£3 <^3(0 = 0. 

To demonstrate that ^3 is integrable and normalized: 

/oo poo r poo 

-00 J —00 L^— 00 



d(3 



let us consider the following affirmation in R : if ai and 02 are two positive real numbers; 
gi is a real function defined and continuous between — ai and ai , and such that for every 
h between — a2 and 02 : 

/ai rai 
gi{x-b)dx= / gi{x)dx; 
-ai J — a\ 

and g2 is a real function defined and continuous between —02 and 02 ; then: 



ai 



ai 



a.2 



gi{x-b)g2{b) db 



a2 



dx = 



ai 



gi{x)dx 



£12 



g2{b)db 



-02 



This affirmation can easily be proved in R by Fubini's Theorem, and so extended to R 
by the VET: 
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If Q!i and a2 are two positive virtual numbers; 0i is a virtual function defined and 
continuous between — cki and a\ , and such that for every /? between —0:2 and 0:2 : 



and 02 is a virtual function defined and continuous between — 0:2 and 0:2 ; then: 



ai 



a2 



.(e-/5)02(/?)ci/? 



a2 



ai 



01 (Oc?^ 



Now, making 



-00 J —00 \.J —82 



dp 



02 



a2 



the above lemma guarantees that taking cki = 00 , 0:2 = £2 , 0i = (^1 , and (f)2 = S2 , we 
get: 

/OO r poo "I r /•£2 

00 L"' — 00 J L"' — £2 

To calculate the integral 



1. 



/OO 
5i(e-/3)<^2(/3-«)ci/? 
-00 

we change the variables n = (3 — a: 

/00—a 
5i{^ - a - fj,)52{iJ,) d/i 
-00—a 



5i{^- a- iJ,)62{iJ,) d/j, 

= 53(e-«). I 

This result shows that the contraction of two continuous Dirac's functions is equivalent 
to a third Dirac's function: 



6i{x - (3)52{(5 -a)d(3 = 5i{x - a). 



VII. Differentiable Dirac's Functions 

We will now consider a differentiable Dirac's function (derivable with continuous 
derivative) and demonstrate the sifting property associated to its derivative. We will 
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say that a real function is differentiable around x E H when there exists an open real 
interval containing x in which the function is defined and differentiable. 

if 5 : R — > R is a differentiable Dirac's function and /: R — > R a real function 
differentiable around a e R, then: 



f 



5'{x-a)f{x)dx= -f'{a). 

Proof: Since 5 is a Dirac's function, there exists a positive infinitesimal e such that: 

\^\>e ^ 5(0 = 5' (0=0. 
By the additivity with respect to the virtual integration interval, we have: 

/oo pa—e /■(!+£ pcxD 

5'(^-a)/(0 = / 5'(e-a)/(0 / S'{^-a)f{0 / S'{^-a)f{^) d^. 
-oo J —oo Ja—s Ja+s 

On the right-hand side of this equation, the integrands of the first and third integrals 
vanish, for any function / , and therefore those integrals exist and are equal to zero. In 
the second integral (right-hand side) we have the product of two continuous functions, so 
the integral on the left-hand side exists and: 

/oo ra+s 
S'{^-a)f{i)di^ / 5'{i-a)f{i)d^. 
-oo Ja—s 

The VET shows that the integration by parts formula holds for virtual integration. 
Since the function /: R — > R is differentiable between a — e e a-|-£,we have: 

S'{^ - a)/(0 = m - a)fmtl - - a)/'(0 d( 

Ja—s Ja—s 

ra+s 

= / S{C-a)[~f'mdC 

J a—s 

/oo 
Sic -a) [-f'mdC 
-oo 

~-/'(a). ■ 

Taking a = and f{x) constantly equal to 1 , we get: 

/oo 
6\x) dx = 0, 
-oo 

therefore: 
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if 5: R — > R is a differentiable Dirac's function then its derivative S': R — > R 
cannot be a Dirac 's function. 

Proceeding on the same line of argument, it is not difficult to obtain the sifting prop- 
erty associated to the higher order derivatives of a sufficiently differentiable Dirac's func- 
tion: 

// 5 : R ^ R is a Dirac's function n times differentiable, and /: R — > R is a real 

function n times differentiable around a G R , then: 

r-oo 



/oo 
5(^)(a; - a)f{x) dx = (-l)^^"^^)- 
-oo 



These sifting properties of the higher order derivatives of a Dirac's function only hold 
for sufficiently differentiable real functions, whereas the Dirac's equivalence, as defined in 
Sec. IV, requires that 

/oo poo 
(f){x) f (x) dx — / ip{x)f{x)dx 
- oo J — oo 

for every continuous real function /: R R . Since there are continuous non-derivable 
real functions, we cannot deduce equivalence formulae in this "strong" sense using those 
properties, even for a Dirac's function sufficiently differentiable. 

Nevertheless, if we are only dealing with differentiable real functions, we can consider 
"weak" Dirac's equivalences, which require integrals equal for this kind of real functions 
only: 

We will say that two integrable virtual functions 0: R — > R and il^-.'K^'R, are 
equivalent in order n when, for every real function n times differentiable /: R — > R , we 
have: 

/oo poo 
(f){x)f{x) dx = tl){x)f{x) dx. 
-oo J —oo 

We will indicate that two virtual functions are equivalent in order n by writing: 

(f){x) = iIj{x) (order n). 

According to this definition, it is not difficult to deduce, from the sifting properties 
above, that: 

// 5 : R ^ R is a Dirac's function n times differentiable, and 51: R — >• R is a real 
function n times differentiable around a e R , then: 

g{x)6^''\x- a) = J2(-'^y(^^9^''~'^(^)^^'^(^-"') (^^^^^ 
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